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The beam coupling impedances of small axisymmetric ob- 
stacles having a semi-elliptical cross section along the beam 
in the vacuum chamber of an accelerator are calculated at 
frequencies for which the wavelength is large compared to a 
typical size of the obstacle. Analytical results are obtained 
for both the irises and the cavities with such a shape which 
allow simple estimates of their broad-band impedances. 



J z (x,y,z;k) = I S(x)5(y) exp(-jkz), 



(2.2) 



I. INTRODUCTION 

High currents in modern accelerators and colliders 
severely restrict the allowed coupling impedance of the 
machine. For this reason, it is important to know the 
impedance contributions even from small discontinuities 
of the vacuum chamber. 

In a recent paper [Q, Kurennoy has analytically cal- 
culated the low-frequency coupling impedance of small 
obstacles protruding into a beam pipe. In this paper we 
present an alternative derivation for an azimuthally sym- 
metric semi-elliptical object protruding into a beam pipe, 
which confirms the dependence on the depth, but not on 
the width, of the protrusion. We also study the more dif- 
ficult — from the analytical point of view — case of an 
axisymmetric semi-elliptical protrusion outside the beam 
pipe (cavity), and present variational results for different 
elliptical eccentricities. 



II. GENERAL ANALYSIS 



Consider a beam pipe of radius i?pi pc and an az- 
imuthally symmetric obstacle whose dimensions are small 
compared with both i? p i pc and A, the rf wavelength. We 
start with the definition of the longitudinal impedance as 



Z [{ (k) 



dvE ■ J*, 



(2.1) 



where the current in the frequency domain for an ultra- 
relativistic point charge is 



*Work supported by the U.S. Department of Energy 
f Present address: AOT-1, LANL, Los Alamos, NM 87545 



with k — w jc = 2tt/5, and with the implied time de- 
pendence of all quantities being exp(jwt). We then iden- 
tify two configurations: the subscript 1 denotes the pipe 
without the obstacle and the subscript 2 denotes the pipe 
with the obstacle. By forming the combination 



J dv(E 2 • J* + x J) 



and using Maxwell's equations to write J in terms of the 
fields E and H, we write the contribution of the obstacle 
to the impedance as 



\h\ 2 Z u {k) 



dS 2 n 2 ■ E\ x if 2 , 



(2.3) 



where the surface integral is only over the surface of the 
obstacle. Using 



™ ZqIq , . 
E lr = - — exp(jfez) 
Zirr 

and 

n 2 dS 2 = 2irr[hdz — zdr], 
with r and z being unit vectors, we have 



Z\i(k) 
Z 



- / drH 2 ^ kz . 



(2.4) 



(2.5) 



(2.6) 



Figure 1 show the geometry for an obstacle protruding 
into and outside of the beam pipe. 
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FIG. 1. The beam pipe with an interior (a) and exterior 
(b) obstacle. One assumes o,()<y,Z<A, -R P i pc . 
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A more explicit form for Eq. (2.6) is 



Zo 



1 



drH 24> {r,z a {r)y kz{ ^ 



Rpipe 



drH 2<i> (r,z b (r))^ kz ^ 



(2.7) 



We now convert the bracket in Eq. (2.7) to a double 
integral for the obstacle in Fig. la: 



-z 



dz — (H 2( p 
oz 



J kz \ 



dr{H 2 ^ kz ) 



—z 



pipe a 

dr / dz — (H 2 ^ zk ) 

Jz b (r) OZ 



-Rpipc 



dr(H 2tj) . 



(2- 



Here Z is a distance large compared with the dimensions 
of the obstacle, but small compared with A and i? p i pe , so 
that H 2t j> exp(jfcz) takes on its value in a pipe without an 
obstacle at z = ±Z, namely 



JT 2 *(±Z) = J exp(TjfcZ). 



(2.9) 



This causes the 2nd and 4th terms on the right side of 
Eq. (2.8) to cancel. We also add the vanishing term 



/ dr / dz—(H 2 ^ kz ), 

Jr p ^-y J-z dz 



and finally obtain 



Z\\{k) 
Zo 



1 

To 



solid drdz—{H 2 ^ kz ), 

border ®^ 



(2.10) 



(2.11) 



where the area of integration is within the solid border 
in Fig. la. Parallel arguments lead to the same result for 
the obstacle in Fig. lb. 

We now use Maxwell's equations to rewrite Eq. (2.11) 

as 



/|(/> ' i - -J!L f L ld drdz(E r - Z Q H^ kz (2.12) 

^0*0 J Jbordor 



Zo 



where we have dropped the subscript 2. For low fre- 
quency we set exp(jfcz) = 1 and obtain 



Z\\(k) 
Zo 



ZqIo 



solid drdz(E r - ZoH^). (2.13) 

border 



Clearly our derivation has resulted in a separation into 
a term involving the electric polarizability and a term 
involving the magnetic susceptibility (in the azimuthal 



direction) as in previous work |3j-|6J. Since the obstacles 
are azimuthally symmetric, we can replace ZoH^ in the 
vicinity of the obstacle by 



ZqH^ — E Q — Z /o/(27ri?pipc) 

and find 

Z\\(k) jk 



Zo 2-KRpipc J VbOTdc 



drdz 



E r — Eq 



(2.14) 



(2.15) 



III. SEMI-ELLIPTICAL INTERIOR OBSTACLE 
— IRIS 

We now proceed to calculate Eq. (2.15) explicitly for 
the geometry of Fig. la. We change from the variable r 
to the variable 



y = R 



pipe 



(3.1) 



with |wl <§; A, \y\ <C i? P i pe , and we use elliptical coordi- 
nates M defined by 



y = c cosh u cos v, z = c sinh u sin v 

a = c sinh Uq , b = c cosh Uq, c 2 = b 2 — a 2 

The metric (Jacobian) is defined by 

dzdy = c 2 D(u,v)dudv, 

where 

D (u, v) = cosh 2 u sin 2 v + sinh 2 u cos 2 v, 
and the Laplacian operator can be written as 

i fd 2 o 2 



d 2 d 2 

+ 



dx 2 dy 2 c 2 D(u,v) \du 2 dv 2 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



The solution to Laplace's equation for the electrostatic 
potential in the region u > Uq, with ty(u ,v) = 0, and 
with the asymptotic field Eq, is 

v) = cE Q cos v [cosh u - e u °~ u cosh u ] (3.6) 

where 

9* 9* 



E r — — 



dr dy 



d 



Eq - cE e u ° cosh u — (cos we u ). 

dy 



(3.7) 



Here we choose cosi> in the second term to preserve the 
symmetry around v = (z = 0), and to satisfy the 
boundary condition for all v at u = uo and for v — ±7r/2 
with u > uq. Also we have 



d du d dv d 
dy dy du dy dv 
1 



cD(u, v) 



sinh u cos v 



du 



cosh u sin v 



dv 



(3.8) 
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Applying this to Eq. (3.7) for E r — Eq, we find from 
Eq. (3.6) that 



gjKjO _ -Jkc 2 



e u ° cosh u Q x 



pipe 

2 . 



dudw [cos 2 ve u sinh u — sin 2 ve u cosh u] 
jkc 2 e u ° coshuo 



2i? D 



soUd duefo [e~ 2u - cos2«], (3.9) 

border 



where we have used Eq. (2.14). 

We now let Z — ► oo and cut off the integration over v 
where 



y = c cosh M max cos t> = F, or 
Y 

Umax = COSh~ > 1 . 

CCOSD 



(3.10) 



This leads to 
z \\( k ) _ jkc 2 e uo coshwo 

Z 47T_Rpi pc 

^ jkc 2 e u ° coshuo 



/ du\e 2u — cosv] 



no 



4tt.R 



pipe 



du cos 2v In 



3 — 2«o 



2F 



CCOS V 



u 



Z\\(k) jkc 2 cosh 2 u a jkb 2 
Z ~ AR ~ 



pipe 



4i? 



(3.11) 



(3.12) 



pipe 



where the last integral in (3.11) was done by parts for 
Y > c. 

For a > b we need to modify our elliptical coordinates 
so that 



y = csinhitsinu, z = c cosh « cos v 

a = ccoshuo , b = csinhwo, c 2 = a 2 — b 2 

The matrix is unchanged, but now 



(3.13) 







1 



dy D(u,v) 



cosh u sin v d sinh u cosv d 



dv 



(3.14) 



This time one finds 

Z\\(k) jkc 2 sinh 2 u jkb 2 



Z 



4-Rpipc 



4i? D 



that is the result is unchanged from Eq. (3.12), again de- 
pending only on the depth of the elliptical protrusion into 
the pipe and not on its width, as also found by Kurennoy 



IV. SEMI-ELLIPTICAL EXTERIOR OBSTACLE 
— CAVITY 

A. Analytical Approach 

We now turn to the exterior semi-elliptical obstacle in 
Fig. lb. Here we need to choose appropriate potential 
forms for u < uq, and for u > uq, and match them at 

U = Uq. 

We again start with Eq. (2.15) and work with the com- 
plete set of solutions of the Laplace equation, namely 
cos nv exp(±nu) and sinra;exp(±nu). For u > uq, (and 
b > a) with the asymptotic field Eq, we choose 



^(u, v) = E y - cE ^2 a n cos ) 



(4.1) 



in order to satisfy the even symmetry about v = (z = 0) 
and *&(u, ±7r/2) = for u > uq. If we write 



#(u ,u) = cE f(v), 



(4.2) 



where f(v) is, as yet, an unknown function, we can solve 
for a n in terms of f(v) to obtain 



\fn + e u ° coshu (5„i, 



where 



/„ = / dvcosnvf(v). 
For u < uq, we write 

oo 

fy(u, v) — (3 n cos tlv cosh ? 



(4.3) 



(4.4) 



(4.5) 



m — 1 



for a potential which is well behaved within the ellipse. 
Recognizing in this case that f(v) = for tt/2 < \v\ < it, 
we solve for (3 m in terms of f(v) to obtain 



A, 



I ri 



7r cosh muo ' 



(4.6) 



where f m is consistent with the definition in Eq. (4.4). 
Both odd and even values of m must be included. 

We now calculate the impedance as we did in the pre- 
vious section, this time including the regions u > uq, 
\v\ < 7r/2 and u < uq, \v\ < tt. For u > uq we find 



(3.15) zl>\k) 



jkc 2 



Z 



47ri? D 



- / L>„ n dudv na n 

>e J Jy<Y n=1 



e -(n-l)u cos ( n + y v _ e -(n+l)u CQs(jl _ 



(4.7) 



Clearly only the terms with n — 1 survive, leading 
ultimately to 
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Z^ik) jkc 2 



-a.\e u ° coshu . 



(4.8) 



For u < uq wc separate the two terms in Eq. (2.15). 
The second is simply 



Z^J(k) jkab j kc 2 cosh u sinh 



U 



Zo 

The first is 



ttAR 



pipe 



AR 



(4.9) 



pipe 



Zf<\k) 



jkc 2 



4:TTR 



dudv 



pipe J Ju<u 



x TO/3 m [cosh(m + 1)ucos(to — 1)?; 

m— 1 

+ cosh(m — 1)mcos(to + l)v]. 
Again, only the term m = 1 survives, and is 



(4.10) 



-/3i cosh m sinh u - (4-11) 



Zo 4i? pipo 
Using Eqs. (4.3) and (4.6) we have for the impedance 



z \\(k) jk 



2i? D 



b ^-+ab-{a + b) 2f —e 



(4.12) 



In order to find /i, we must obtain and solve the in- 
tegral equation which represents the match of d^S / du at 
u = uq, \v\ < 7r/2. Here 



du 



cE 



U — U o + 



sinh uq cos v 



na„e ™"° coshu 



(4.13) 



and 



9zt 



c_Eq w/3 m sinh miip cos mv. (4.14) 



m— 1 



Equating Eqs. (4.13) and (4.14), and using Eqs. (4.3) and 
(4.6), wc find 

j 2 dv'f(v')K(v,v') = ircosve u ° (4.15) 

where 

CO 

/T(w,t/) = y^(2 + tanhnuo) jicostocos ™' 

?i=i 

odd 
CO 

+ m tanh roitp cos mv cos m«. (4-16) 



We now multiply Eq. (4.15) by dvf(v) to obtain 



dvf(v) cosv 



n2 



Jj„ dvf(v) cosv 



/J } *>/(") /J } dv>f(v')K(v, v') Ke u ° 

fie- ua 



(4.17) 



This is a variational form for /i, the only unknown pa- 
rameter in Eq. (4.12) for the impedance. An accurate 
numerical value for Z\\ (k)/Z can be found by expanding 
f(v) into a complete set in the interval \v\ < 7r/2, then 
truncating and solving the resulting matrix equations ob- 
tained by maximizing Eq. (4.17). We write 



/(«) = Xl c P sin 



7T\ COSpV 



P =l 

odd 



p 



(4.18) 



truncated at p — P and normalize f(v) so that C\ = 1. 
This leads to 



/ie-"° 



-ffn + 2 ^ c p H p i 



P=3 
odd 



p p 



p — 3 q — 3 
odd odd 



(4.19) 



where the symmetric matrix i? pg is 
2 + tanhpwo r 



16 x 



m tanh muo 



7T 2 ^ (m 2 — p 2 )(m 2 — q 2 ) 



(4.20) 



Maximizing Eq. (4.19) with respect to the coefficients c p , 
p = 3, 5, • • • , P, leads to 



fie- 



TT 



-i -1 



P P 



Hn - ^2^2 Hi p( h 1 )pi h qi 



p = 3 q = 3 
odd odd 



• (4.21) 



Here (H 1 ) pq is the inverse of the matrix H pq with p and 
q = 3, 5, 7 • • • P. This square matrix has the dimension 



(P-l)/2by (P-l)/2. 



Note that 



and 



tanhpwo = (1 - w p ) / (1 + w p ) 
tanhmwo = (1 - w m ) / (1 + w m ) 



(4.22) 
(4.23) 
(4.24) 
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with 



w=(b- a)/(b + a) 



(4.25) 



The final result for the impedance is given in Eq. (4.12), 
using Eqs. (4.20) and (4.21). 

The analysis for an obstacle with a > b proceeds in 
a similar, but not identical pattern. The result is once 
again Eq. (4.12), using Eqs. (4.20) and (4.21), with only 
one change in Eq. (4.20): 

tanhpito — » cothpuo, 

with tanhwo now being a/b instead of b/a. Thus w in 
Eq. (4.25) is replaced by — w in Eq. (4.23), but the use of 
cothpuo instead of tanhpuo for odd p leaves the expres- 
sion for H pq in terms of a and b unchanged. The same is 
true for the term tanhTOUo m Eq. (4.24) since m is even. 
So the final expression in Eq. (4.12) is unchanged pro- 
vided H pq in Eqs. (4.20) and (4.21) is expressed in terms 
of a and b. 



B. Variational Approach - Numerical Results 

We proceed with a numerical investigation of the vari- 
ational scheme described by Eqs. (4.12)-(4.21). Truncat- 
ing the sum in the denominator of Eq. (4.21) at differ- 
ent N = (P — l)/2 = 1, 2, 3, . . ., we explore the scheme 
convergence, and compare the results for the impedance 
(4.12) with those obtained by other methods. In doing 
so, it is convenient to rewrite Eq. (4.12) in the following 
form 



Z\\ (k) jk nab (a 
-r 



(!) 



where 



1 



Zq 27r_Rpi pc 2 

2(l/a; + 2 + a;) 



[ ' x 2 + x + l&s{x)/n 2 -^{x)' 



and 



(1 + x) 2n - (1 - xf 



^ OO 

^ = ^ S (" 2 - V4) 2 (1 + x? n + (1 - xf 



(4.26) 



(4.27) 



(4.28) 



Here £(x) denotes the sum in the denominator of 
Eq. (4.21) which is to be truncated. 

The advantage of the representation (4.26) is that we 
know the asymptotic behavior of F(x) for two limiting 
cases. For x <C 1, i.e., when a -C b, but still b <^C i? p i pe 
— a short and deep enlargement — it has been demon- 
strated in B that 



F(x) 



1 



(4.29) 



In this limit, the inductive impedance in Eq. (4.26) is 
mostly of magnetic origin: the beam magnetic field fills 



the cavity volume without being substantially perturbed, 
and therefore the inductance is simply proportional to 
the area of the obstacle cross section. A correction of the 
order of x — a/b to this term comes from the electric con- 
tribution. For a deep pillbox of depth h which is much 
larger than width g, the electric contribution was calcu- 
lated in |^| by means of conformal mapping. It results 
in the electric term — g 2 /(2n), which is small compared 
to the magnetic one, equal to hg for such a pillbox. Ob- 
viously, the shape of a short and deep enlargement - 
rectangular or semi-elliptical — does not affect the elec- 
tric term as long as g <C h, since the beam electric field 
does not penetrate deeply into such a cavity, unlike the 
magnetic one. Substituting g = 2a into the electric term, 
and replacing the pillbox area hg by the semi-ellipse area 
irab/2 leads to the asymptotic form in Eq. (4.29). 

The opposite limit, x 3> 1, corresponds to a very shal- 
low cavity, a <C b. It has been shown for many particu- 
lar shapes of such cavities (see Q and references therein) 
that the low-frequency impedance of a small shallow cav- 
ity of the depth h and of an iris with the same cross 
section and having the same depth, are both inductive, 
equal to each other, and in the leading order are pro- 
portional to h 2 . Since we already know the answer for 
a semi-elliptical iris (see Section III) we expect that for 
x > 1 



F(x) 



1 



(4.30) 



to match the low- frequency impedance of the shallow iris, 
given in Eq. (3.15). 




FIG. 2. Function F(x) versus ellipse aspect ratio x = a/b. 
The solid curve is an interpolation of numerical results (thick 
points). The dashed lines show the asymptotic behavior Eqs. 
(4.29) and (4.30). 

The results of our numerical study are shown in Fig. 2, 
where the function F(x) is plotted against the ellipse 
aspect ratio x = a/b. The convergence of the variational 
scheme is rather fast for all values of x; in fact, results 
obtained with N — 1 (i.e., when the matrix in the sum 
E(cc) is truncated to merely a single number) and those 
for N = 8, when the matrix has size 8x8, differ by less 
than 0.5%. And, of course, we can obtain the asymptotic 
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value of F(x) for large N with much better accuracy, 
well below 10 -3 , simply by extrapolating the results for 
different matrix sizes at fixed x. Figure 2 also shows very 
good agreement with the expected asymptotic behavior 
Eq. (4.29) for small x and Eq. (4.30) for large x. 



V. TRANSVERSE COUPLING IMPEDANCE 

We start with a dipole drive current for the transverse 
impedance in the form [pf 



J z (x, y, z; k) = I 5(y) exp(-jfcz) x 

[S(x - xi) - S(x + xi)], 



(5.1) 



where we eventually proceed to the limit xi — » 0. 
It is straightforward to show || that the transverse 
impedance in the x direction can be written as 



Z x (k) = - 



1 



4fcx?|/ | 2 



dvE ■ J* 



(5.2) 



analogous to Eq. (2.1) for the longitudinal impedance. 
Use of Maxwell's equations as we did in Section II, leads 
to 



Z x {k) 



1 



Akx{\I \ 2 Js^ Sl 



dS z n 2 ■ E\ x H 2 , (5.3) 



but we must now use the form of E\ (and H 2 ) appropri- 
ate to the source current in Eq. (4.26). In fact, we now 
have for E\ and ZqH\ at the beam pipe wall 

2Z I 

Ei r = ZqHi^ = — — 2 — x\ cos(/)exp(— jkz) (5-4) 



where Eq, the maximum asymptotic field at the wall, is 
2ZqIq 



En 



(5.8) 



Comparison of Eq. (5.8 ) with Eq. (2.15) shows that 
the calculations for an exterior and an interior obsta- 
cle are exactly the same as they were for the longitu- 
dinal impedance. In fact, the results for the transverse 
impedance can be obtained simply by multiplying the 
results for the longitudinal impedance in Eqs. (3.12), 
(3.15), (4.12) and (4.21) by 2/kR 2 pipe . 
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Eitj, — ZqHij, — 0. 
As a result, we can write 

Z x {k) „ 1 



Zn 2kxiI ZnirR p i pc 

jjcostj) J dr(ZoH 2 4,e 3kz ) 



(5.5) 



Once again we have written the impedance as an integral 
along the surface of the obstacle, where H 2< f, arises from 
the driving field components E\ r and El\$ at the wall. 
Dropping the subscript 2, extracting the factor cos <j> from 
and integrating over <j) leads to 



Z x (k) 
Z 



1 



2kxiIoZ R p i pc 



dr{Z H^ kz ). (5.6) 



We now write Eq. (5.6) as a double integral over drdz 
as we did in Section II, obtaining 



Z x {k) 



-J 



d3 / solid 
pipe -'border 



dr I dz 



E r — Eq 
E 



(5.7) 
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